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5— ( ' We propose a simple phenomenological model describing composite crystals, constructed from 

, ^ I two parallel sets of periodic inter-penetrating chains. In the harmonic approximation and neglecting 

thermal fluctuations we find the eigenmodes of the system. It is shown that at high frequencies there 
are two longitudinal sound modes with standard attenuation, while in the low frequency region there 
is one propagating sound mode and an over-damped phase mode. The crossover between these two 
regions is analyzed numerically and the dynamical structure factor is calculated. It is shown that 
the qualitative features of the experimentally observed spectra can be consistently described by our 
model. 

o 

C/3 , PACS numbers: 61.44.-n , 63.20.Dj 

I. INTRODUCTION 

' Among the various types of incommensurate systems one of the simplest kind (the so-called uniaxial composite 
crystal) is constructed from two sets of regular inter-penetrating incommensurate chains. The recurrent interest for 
these systems is related mainly with the experimental activity on several classes of intergrowth compounds which can 
be viewed as good physical realizations of the uniaxial composite crystal. There is a considerable literature (mainly 
Q • theoretical but not only) discussing eigenmodes and related properties of composite systems (see e.g. the review 
O I article |^ and references therein, and the more recent papers quoted in ||^, |^). 

Although a number of sophisticated calculations have been published over the last 20 years ||, ||, ^ ||, ^ , there 
is still a clear need for a simple (but yet non-trivial) theoretical model with predictions which can be directly tested 
^ 1 experimentally. In previous models all branches of excitations were simultaneously examined and many subtle 
[ details of the composite systems were discussed in connection with higher order commensurability effects Q , such as 
the hierarchical nature of gaps arising at commensurate - incommensurate phase transitions ^ . The results of these 
' investigations, even in the long wavelength limit, are expressed in terms of many unknown parameters, and lead to 
many different modes, while experimentally it is not clear whether it is possible at all to observe these modes and to 
CS| [ determine these parameters. Furthermore, although the basic theory of excitations in incommensurate and composite 
. systems was born long ago , [0 , , many questions of principle remain unsettled, and there is an impressive 
quantity of unexplained, partially explained or contradictory results. In part this frustrating situation is just due to 
the lack of a simple and tractable analytical model. 

Our motivation for adding one more paper to the topic is precisely to propose such a model. We do not intend 
to discuss in the present paper any of the dynamical details connected with the description of discontinuity edges 
P5 I in the vibrational density of states nor do we deal with the many structural properties (see e.g. |^) related to 
O ■ commensurate - incommensurate phase transitions. Our aim is much more modest and based on the fact that many 
robust and experimentally testable features of the excitation spectra are not sensitive to the more delicate aspects of 
commensurability effects. 

Note that one can consider composite structures as a new state of matter like liquid or solid states. Composite 
structures differ from modulated crystals insofar as one can not define an average periodic structure for the full system. 
Unlike crystals with incommensurate density waves that are characterized by weak deviations of atoms relative to 
their regular positions and unlike the adsorption structures with the rigid host and flexible guest sublattices, both 
subsystems in composites are more or less equivalent and therefore one can not treat them perturbatively as an 
external potential acting on its elastic counterpart, nor can one use the classical Frenkel - Kontorova model 
to describe phase diagrams and excitations. Note that composite incommensurate crystals are also different from 
quasicrystals. According to the traditional classification quasicrystals refer to systems where the rotational 
symmetry corresponds to a forbidden crystal symmetry, which determines the unique ratio of incommensurate length 
scales that defines the structure, whereas for incommensurate crystals the rotational symmetry corresponds to an 
allowed crystal symmetry and many different sets of incommensurate length scales are possible. 
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Due to the existence of (at least) two length scales, there are twice as many hydrodynamic variables as in conven- 
tional crystals. However in quasicrystals these additional variables cannot be interpreted as real displacements of any 
of the sublattices, and this distinction leads to specific physical consequences. 

The most interesting aspect of the physics of composites is the dynamics of their sublattices. In order to formulate 
the problem, let us consider first two inter-penetrating atomic chains with periods a and b and assume no coupling 
between them. Then, two independent sets of incommensurate (in g-space) phonon branches with dispersions Wa(<z) = 
Cafl"^ |sinq'a| and ujb(q) = Cbb~^ \sinqb\ should be observed (Fig.^. What happens with them, when the inter-chain 
interaction is properly taken into account? Different types of coupling are responsible for different effects. 
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FIG. 1: Phonon spectra of two uncoupled incommensurate chains a and b. 

The elastic coupling results in the renormalization of the sound velocities |^ . It leads to repulsion of the crossing 
branches with formation of intermediate gaps. The umklapp intersection of branches occurs however for any q since, 
because of the period incommensurability, the condition uja{Nq) « ujb{Mq) can always be satisfied when the integers 
N,M are chosen properly. Hence a rich hierarchical set of gaps in phonon spectra is expected (see e.g. [^). Some 
general features related to breaking of analyticity due to static elastic couplings in the Id double-chain model were 
reported recently in 

The electrostatic coupling takes place when the chains a and b are charged. It leads to nonlocal renormalization of 
the phonon self-energy, nontrivial retardation effects, and to the plasmon gap formation [D, Q. 

The dynamical dissipative coupling originates from the internal friction between the chains. To our knowledge the 
dissipative dynamics of composites was studied only by the authors in Refs. Q and In the papers ^, ^ the 
authors have taken into account dynamic resonance processes arising due to the existence of degeneracy points of the 
subsystem modes which allow the interchange of energy and momentum. They described these processes by including 
into the response matrix a finite (in the limit g — > 0, and w — > 0) off-diagonal contribution. For a strictly harmonic 
lattice the contribution vanishes identically. We propose below a microscopic realization of the dissipative dynamical 
coupling which is always present in real systems even in harmonic approximation. As it concerns to the paper |^ 
the authors took into consideration all possible types of couplings and all degrees of freedom of composite systems, 
and therefore it is difficult to survey their results (and all the more to compare the results with experimental data) 
because of the large number of hydrodynamic variables involved. Meanwhile the broadening and renormalization of 
phonon spectra due to the inter-chain friction can easily hide the subtle hierarchical gap structure as well as other 
minute effects, like the fluctuation - driven destruction of the Id composite long - range order or the static and 
dynamic lock-in effects |2). 

It is the purpose of the present publication to provide a simple and tractable model that incorporates the dissipative 
dynamic coupling in composite systems to calculate the experimentally testable features of the excitation spectra. 
We put forward a model of a composite system, which can be used to find the eigenmode spectrum (under relatively 
weak assumptions), and to compute the dynamical structure factor. 

The outline of this article is as follows. In Section II we formulate a model description of a composite incommensurate 
crystal which includes all the basic dynamical equations. In Sections III and IV the calculations of the eigenmodes 
and dynamical structure factors are presented. Section V is devoted to a discussion and summary of our main results. 
The Appendix deals with a very schematic description of a particular microscopic model of a uniaxial composite 
system. 
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II. MODEL 

The simplest composite crystal consists of alternating chains with incommensurate periods a and b that are located 
in the x — y plane and are directed along axes parallel to x aX y — md and at y = (m + l/2)rf. In this case even the 
basic structure is not commensurate. Because of the local differences in environment the a chain is modulated with 
period 6, and the b chain - with period a. Therefore the atomic positions can be given by 

R'^(n, m; t) = mdey + nae^ + u"(n, to; t) + 
e^gi{na - A) ,{n,m = 0, ±1, ...) ; 
R^(n, m; t) ^ {m + l/2)dey + nbe^ + u^{n, m; t) + A + 

e^g2{mb + A), (1) 

where A is a certain natural offset of one set of chains with respect to its counterpart, VL°-'^{n, to; t) are the displacements 
of the atoms about their average equilibrium positions, gi{x + b) = gi{x) and g2{x + a) = g2{x) are modulation 
functions. 

In fact just the definitions (|l|) allow us to introduce all the degrees of freedom relevant for our model of a composite 
system (including therefore - the main branches of excitations). If we neglect for a moment the interchain interaction 
(purely hypothetically since in the absence of subsystem b subsystem a would be generally unstable), we get evidently 
two zero - energy modes for the two chains, each one implying a rigid displacement of one chain with respect to its 
fixed counterpart. 

Interchain interactions introduce new ingredients in this scheme. In the incommensurate state the position of the 
centre of mass of the atoms relative to the modulation potential can be described by a phase variable. In a rigid 
chain approximation the acoustic phonon mode corresponds to a global rigid displacement of the two sets of chains, 
while the phase mode corresponds to relative rigid displacements of the two chains in opposite directions such that 
the centre of mass of the composite crystal remains fixed. In a commensurate state this phase variable is pinned, 
thus there is one acoustical mode and one optical mode with a gap whose magnitude is determined by the pinning 
potential and strongly decays with the order of commensurability. In the latter case our consideration will be still 
valid for characteristic frequencies higher than the gap. 

Considering the relative displacements to be rigid means that we neglect the above mentioned intermodulation. 
The interchain interaction moves the atoms on both chains away from their regularly distributed positions. However 
if the modulation functions gi and 52 are differentiable, then the atomic displacements corresponding to the phase 
mode are given by the following "dressed" displacements u 

u°(n, to; t) u''(ri, to; t) + e^bg^ , u^(ri, m; t) u''(n, to; t) + exag'2 , (2) 

where g[ denote a derivative over the argument. 

In order to discuss the long-wavelength, low-frequency excitations one should transform the discrete variables intro- 
duced above into continuous ones. For sound modes which consist of equal shifts of R"(ri, m;t) and R''(n,m;t) this 
transformation is evident. Namely in a continuum approximation the discrete variables \i°''^{n,m]t) are substituted 
by continuous functions u°'''(r, t). However it is not the case for the phase mode which consists of inhomogeneous (on 
a microscopic scale) atomic displacements (corresponding to sliding of one chain system with respect to the other). 
However when the modulation functions gi are analytic, one can still use the long- wavelength description but for 
dressed displacements u. 

For the benefit of the more skeptical reader the question concerning " dressed" variables should be clarified in more 
detail. In fact the dressed displacements u depend on the displacements and u** as initially defined according to (|l|) 
(since the modulation functions depend on the difference u° — -u^). For weak enough interactions and far enough away 
from commensurability, the modulation functions gi^2 are continuous analytic functions and the energy of the system 
remains unchanged when the phase m° — is varied. However a uniform shift — does not indicate a uniform 
displacement of the atoms in both chains, since the functions gi and 172 are not constants. Thus strictly speaking 
we should write down equations of motion separately for dressed and undressed variables. In our phenomenological 
approach it means a certain renormalization of all phenomenological constants entering the equations. In what 
follows we neglect the renormalization and will use u°'^(r, t) for both (acoustical and phase) modes. In spite of 
this erroneous assumption (which is equivalent to 51^2 — const) the approximation correctly identifies the important 
modes and characteristic scales in the problem. This sin of omission can be easily corrected when detailed information 
concerning the values of the parameters becomes available. 
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The elastic energy can be written as a quadratic form of the strain tensor efj = ^{diu" + djuf) : 

^ = \K^H<A + r^<~<f-^ (3) 
a,/3 = a, 6; i,i,k,l = x,y. 

The inter-chain interaction is given by the relative displacements along y (the 2nd term in (^) and also by the cross 
terms in the elasticity tensor Xfj/^i- At the same time the free energy expansion is invariant with respect to the 
relative shift of the chains along x because of the incommensurability of the two periods, therefore there is no terms 
depending explicitely on u"''' (only on gradients of the displacements, i.e. on e"'''). 

As a note of caution we should remark here that there is a fundamental difference between the both kind of soft 
degrees of freedom in (|3|), i.e. between acoustic and phase modes. The first one is a true Goldstone mode, because it 
breaks the continuous translational symmetry (and therefore remains hydrodynamic and gapless even if the modulation 
functions are not continuous). But the latter is a pseudo-Goldstone mode and it loses its hydrodynamic character in 
the case of discontinuous modulation functions (since phase variable breaks only countable translational symmetry). 

To complete the derivation of the dynamical equations we have to know not only the elastic energy (^ but also the 
dissipative function R. The internal friction produced by the sliding displacement of the chains results in a dynamical 
inter-chain coupling proportional to the chain velocity difference — Uj. that contributes to the dissipation function: 

D ^rr <^ ■b\2 I -'- Q/3 -a -0 /.N 

2 i^x -uj + ^Vijki'^ij^ki ■ (4) 

The second term in (^ corresponds to the viscous dissipation, and the coefficients ry^^; entering (^ have the standard 
meaning of viscosity tensor elements [p^ . One new phenomenological coefficient F describes a friction proportional to 
the velocity difference between the two chains. A few comments to assign physical meanings to the terms in equation 
(^ may be helpful to the reader. The first term in (^ is the energy dissipation caused by relative motion between the 
two set of chains. In the framework of our purely phenomenological approach we are not in a position to determine 
a dominant microscopic mechanism for the dissipation. A natural estimate for the order of magnitude of the friction 
coefficient F (based on the assumption that the work produced by the drag force density is completely dissipated in a 
viscous media) is F oc rj/dP', where 77 is the characteristic viscosity. It is worth noting that dynamic dissipative, static 
elastic and lock-in interchain interactions can differ both quantitatively, since there is no simple relation between the 
coupling coefficients F, cr, and the elastic moduli, and qualitatively, because unlike static interactions, the dynamic 
dissipative coupling does not necessarily produces intermodulations of the two subsystems. 
Now we are in a position to write down the equations of motion 

..a dT OR 



where 



Pa/5 = ( ) (6) 

is the matrix of densities. Eq. (^ is the second Newton law where the forces acting on the component a are: (i) the 
friction between the components due to their relative motion, and (ii) the elastic forces. 

As it is well known the propagation of elastic waves in anisotropic media is a rather complicated phenomenon. 
The directions of polarization and the eigen-frequencies are determined by the dispersion equation which is derived 
for the monochromatic solutions of (|^). In a general anisotropic case however none of the polarization directions 
correspond either purely longitudinal or purely transverse directions with respect to the direction of the wave vector. 
We can look for solutions of the solution of the equations of motion (|^) in the form 



a i(qr-wt) 



(7) 



The dispersion ijj{c\) can be found by substitution of ^ into (|^) and, then, by diagonalization of the corresponding 
dynamical equations: 

^a'/3(9,C.)<o = 0. (8) 
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The four degrees of freedom of uf {a = a,b, i = x, y) correspond to the four propagating eigenmodes. 

We are interesting in the case when the excited wave is propagating along x (other geometries wih be shortly 
discussed in the Conclusion). Then, the equation of motion (||) can be split into a transverse part for Uy and a 
longitudinal part for u". In the transverse case, the last term in Eq. (^) results in the decoupling of the excitations 
into sound and optical modes. Thus for the polarization perpendicular to the chains we will find the usual optical and 
acoustical branches. Less trivial is the longitudinal case when the last term in Eq. (|^) vanishes and the dynamical 
mode coupling T5ix(ui —Ui) provides the most important contribution to the properties of the propagating excitations. 
In the next Section we consider this case in detail. 



III. LONGITUDINAL EXCITATION SPECTRUM 



The dynamical equation for the longitudinal eigenmode with q || x looks like 
with 



(9) 



\ab 2 

' ,7. C/ 



iTuj 



iVuj - i-nfxxxT^ 



,2 I \hb 



Kxxx<i' + 



iVxxxxQ^^ 

2 _ ,„bb „2 



q uj 



(10) 



The terms with F and A"'' provide the dynamical dissipative and clastic coupling between two sound waves prop- 
agating in chains a and b. To pick out the effect of this coupling on the phonon dispersions we neglect the intrinsic 
phonon attenuation that is given by terms with rj. 

To find the eigenmodes one should solve the characteristic equation: 



detDmq,uj) = Pc,Pb[{uJ^ - clq^ + ilaUj)iuJ^ ~ clq^ 
+i76w) - {cfq^ +i{'yalby^^ujf]=0 



(11) 



where: 



2 _ \aa I 2 _ \bb 1 2 \ab 11 n1/2 

~ ^xxxxl Hon <^b ~ ^xxxxl Pb, i^i — ^xxxxl \HaPb) , 

la^T/Pa, Ib^^/Pb, 



(12) 



At small q (q < ^a.bl(^a,b) the dynamical dissipative coupling substantially renormalizes the dispersions of the two 
sound modes a and h. 

The first two branches correspond to the time-conjugated sound-like modes 



with effective velocity and damping: 



1 



-facl+-fbcl-2{jajb)^^^ 
7a + 76 

{la-lb)c^~ilalbf' (cl-cl) 



(13) 



(14) 



Vs 



2 (7a + IbV laCl + Ibci - 2 (7a7;>)'/' cj 

The third mode corresponds to the over-damped phason-likc mode with pure imaginary dispersion (diffusive mode): 

UJ « -iTJphq'^, (15) 

where the damping factor is given by: 



.2 



a'-b 



laCl + IbCl - 2 (7a 76) 



1/2 J2 ■ 



(16) 
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In addition there is a purely relaxational fourth mode: 

Ul = -i{ja + lb) (17) 

which, however, is hidden by the third mode that relaxes much more slowly. 

At high q: {q > ia,b/ca.b) the elastic coupling leads to the effective renormalization of the velocities of the sound 
modes a and b. The effect of the dynamical (dissipative) coupling is manifested by a specific q -independent sound 
attenuation. The dispersion of these two time-conjugated sound modes reads: 

uj = ±ci2(j — 1712- (18) 

The renormalized velocities and attenuation factor are given by: 



2cl,^cl{l±K) + cl{lT^), (19) 

4712 = [Vla{l ± K-1) - v/7fc(l T K-1))^ . (20) 
where the elastic coupling parameter k: 

is larger than one (weak coupling) and smaller than (c^ -I- c^)/ — (stability condition). 

At higher q the intrinsic phonon attenuation ^ ^"iiVxxxx/ Pia,b))(f is superimposed onto the dynamical attenuation 
—1712, a situation which can lead to a nontrivial (/-dependence of the phonon line- width. 

Of course it is possible to find the eigenmodes (at least no problem at all to find the dispersion laws numerically) 
in the general case of arbitrary q taking into account all phenomenological parameters entering (|lo|). However aiming 
for a simple model we have chosen the simplest coupling form for the displacement components on both chains, 
following the principle of minimal requirements. To keep our model tractable (but yet non-trivial) and to avoid 
tedious calculations we drop the intra-mode coupling X'^^xxQ^ that does not change the qualitative picture. As will 
be discussed further, the intrinsic phonon attenuation expressed by iri'^^,j.^q^u) can be also neglected in the low energy 
region which is relevant to this discussion. Although actually this picture is not entirely correct it correctly identifies 
the important modes and characteristic scales in the problem. 

As a note of caution we should also keep in mind that there are four types of interactions involved, namely: 

• the elastic interaction proportional to gradients of the displacements, and related to off-diagonal terms of the 
elasticity tensor; 

• locking interaction, proportional to sublattice displacement differences (the last term in (^) with the coefficient 

• two dissipative couplings described by the off - diagonal terms of the viscosity tensor and by the mutual friction 
coefficient F in (||). 

The most robust phenomena are related to the dissipative couplings which we have taken into account in the 
derivation of the dispersion laws (|l3|- ^8|). However the locking and elastic couplings could lead to appreciable effects, 
like e.g. a gap at g = for the sliding mode. 



IV. DYNAMICAL STRUCTURE FACTOR: SCATTERING EXPERIMENTS 

The above discussed distinctive features of the excitation spectrum should be observable in different kinds of scat- 
tering experiments: ultrasonic measurements, Brillouin light scattering and inelastic coherent scattering of neutrons 
or X-rays. The double - differential cross-section gives the probability that say a neutron with initial energy E is 
scattered (coherently or incoherently) into a detector subtending a solid angle dTl around a certain direction with 
spherical coordinates 29 and </>, with a final energy E' . We note hu) = E — E' , {Tiuj is the energy transfer), and 
?iQ = ?ik — Tik! is the momentum transfer, where ?ik and Tik! , are the incoming and outgoing neutron momenta. The 
scattered intensity is proportional to the differential cross section for inelastic scattering: 
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dVLdio k 2Trh 

where angular brackets denote a thermal average, Sa = Sa, Sb are the coherent atomic scattering factors i.e. differential 
cross sections of light, sound, neutrons or X-rays for single atoms of chains a and b, and the integrand in (^2|) is called 
the dynamical structure factor ^(Q,^) which is accessible by scattering techniques, and obeys the so-called detailed 
balance condition 



5(Q,c.) = exp( ^)^(-Q,-c.) 



if the sample is in thermal equilibrium. 



According to the fluctuation dissipation theorem |16 



1 — exp(— aw/ J j TT 

where x" is the imaginary part of the generalized susceptibility. 

Calculations of x" is a formidable task that can be done only on the basis of an appropriate microscopic model 
(see e.g. the Appendix to our paper). In a general case the dynamic susceptibility x is a matrix with respect to 
the displacement fields of the two chains and to the respective Cartesian components. Owing to the lack of discrete 
translational invariance in the x - direction it is non-diagonal also in the wave vectors. Since the full investigation of 2d 
(all the more 3d) problem requires an extraordinary amounts of computational work we must defer the investigation 
of the full problem to the time when experimental results suitable for a quantitative comparison become available 
(see also the Conclusion section of the paper), and present here only the results for the wave vector parallel to the 
incommensurate direction {x - axis). 

In the case of longitudinal phonons it can be written in the following form (we dropped here unessential factors, 
since we are not interested in the absolute value of S'(Q,w) or x"(Q,'i^)): 

S{Q,u)o,Re J2 SaS;3e'«("""''"')Q'^/mxS^(n,n',c^). (23) 

n,n' ,a.J3 

To study the structure factor we need to know the equilibrium positions of the atoms, the amplitude of oscillations, 
and the dispersion relation. Thus the generalized susceptibility Xa/ji which is the response function, is proportional 
to the inverse dynamical matrix 

Xr^(n,n',w)^DS^(n,n',w)-i. (24) 

We discuss here the cases when the momentum transfer Q (parallel to the x -axis) lies close to the Bragg peaks of 
the composite crystal Q — Q™^" + q with 

where m, n are integer numbers. Expression ( |25| ) implies that each observed Bragg peak can be classified as follows: 
main reflections of the a subsystem (m 7^ 0, n = 0); main reflections of the b subsystem (to = 0, n 7^ 0); common 
reflections to a and b subsystems (to = 0, n = 0); pure satellite reflections (to 7^ 0, rt 7^ 0). With regard to the general 
program we do not have the ambition to find the dynamical structure factor in the whole region of parameters. In 
fact such a general formulation would not be tractable and in any case it would have little meaning in view of the 
qualitative nature of our model. For this reason and only for illustration purposes we consider a few particular cases. 

Consider first the forward scattering when the momentum transfer Q = q is small (acoustic case) . This regime cor- 
responds to Brillouin (light and X-ray) scattering. On the basis of ( p3| ) the structure factor S{Q, ui) is reduced to (aim- 
ing at a qualitative description of the dynamics we drop the amplitude factor (fc'/7rfc) exp[— 2IVq] (1 — exp(?icj/T)) , 
where Wq is the Debye - Waller factor): 

S{q,Lo) cx J2ilV^Wsi3lmDl^fiiq,cu)-' (26) 

a/3 
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FIG. 2: The structure factor of a uniaxial composite system. Bold lines correspond to S{q,ui), normal lines correspond to 
S(Q^ + q,ij-i) and S{Q]!^ + q,ui)- Ticks: "[" and "]" correspond to Caq and c^^q; ticks "|" and "x" correspond to c+g and c_g. 
Parameters: pa = 1, pi, = 2, Ca = I, Cb = 2, 7a + 75 = 0.45. 



where D^p{q,uj) is given by (|To|). 

As mentioned earlier in order to pick out the effects of the dynamical dissipative coupUng we neglect the elastic 
coupling A"'' and assume that both channels a and b are equally-participating in the scattering process { Sa = si, — s). 
Then S{q, u) can be written in a simple form: 

S{q,uj)(X (27) 

5^ (7a + lb){Pa - Pbfq^[<Jj'^ - Clq^? 

(w2 - c2g2)2(t^2 _ c2g2)2 _|_ _,_ -y^)t^2(-^^2 _ c2_g2)2 ' 

with: 

c± — T ■ 

Pa ± Pb 

The structure factor in the so-called pseudo acoustic regime (when Q is close to a nonzero Bragg peak) is different 
from the acoustic structure factor (p7|). It is worth noting also that this feature is different from the case of phonons 
in conventional crystals. 

Consider e.g. the pseudo acoustic region close to the Bragg peak Qia = 7r/a (Fig. 1). Because of the large energy 
separation, the modes a and b are decoupled and the structure factor for mode a is written as: 

Sa{n/a + g,c.) cx {Q,Ju;)slD:^{n/a + q, cj)-\ (29) 

The matrix element D^^in / a + q, lu)^^ coincides with D^^{q, u;)^^ (it corresponds to the umklapp coupling processes). 
On one side the "a" phonon with Q = n/a + q is decoupled from the "6" phonon with Q = ir/a + q because of the 
large energy separation. At the same time the "a" phonon with Q = n/a + q is equivalent to the "a" phonon with q, 
that in turn, is strongly coupled with the "6" phonon with Q — q. 
A straightforward calculation gives 

5a(^+9,Cc>) cx (30) 

QUla+li,)pl{u^'~cWf 

(w2 - C2g2)2(^2 _ c2g2)2 _^ (7^ + 7b)cj2 (^^2 _ c\q^Y ' 



The results of the calculations of ^(q, ui) for different values of u! scanning the structure factor at a fixed value of q 
are shown in Fig. 2. The picture presented above leads to the following qualitative predictions: 
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• The maximum in ^(q, w) due to the sublattice phonon with the lower frequency disappears at a certain wavevec- 
tor. 

• Both S{q,uj) and S{m:/a + S(rmT/b + q,Lo) have asymmetric profiles with " antiresonance" points where 
they vanish (as in the case of the Fano resonance, when the interference of a discrete state with a continuum 
gives rise to characteristically asymmetric peaks in excitation spectra); 

• Although S{m:/a+q, lu), S{'rm:/b+q, uj) correspond to "decoupled" phonons. they show traces of mixing between 
the two phonon branches because of the umklapp process; 

• At low q we always have a double peak structure with a central peak and a sound-like side peak. 

V. CONCLUSION 

Let us sum up the results of our paper. We propose a simple phenomenological model describing a composite crystal, 
constructed from two sets of periodic inter-penetrating chains having a common direction (cither commensurate or 
incommensurate). We intentionally restrict ourselves to the simplest case taking into account mainly the dissipative 
coupling of the two subsystems, since theoretical modelling of composite crystals is often hampered by lack of detailed 
knowledge concerning the parameters entering more complete and more sophisticated models than the present one. 

As incommensurate composites do not possess translational symmetry, one cannot define a Brillouin zone. In zero 
order approximation if both sublattices are independent, there are two phonon modes (one for each sublattice) but only 
one sound velocity (in the long- wavelength limit). Coupling between sublattices leads (again in the hydrodynamic, 
i.e. q limit) to one acoustic mode for the composite, and one sliding mode, as found above. In this limit the 
acoustic mode is associated with a global translation of the composite, while the sliding mode describes the relative 
"dressed" translation of the two sublattices. 

We find the eigenmodes of the system and calculate the dynamical structure factor S{Q,uj). We show that at 
high frequencies there are two longitudinal sound modes with wave vector independent attenuation, while in the 
low frequency region there is one propagating sound mode and two over-damped phase modes. The structure factor 
has a specific profile resembling the Fano resonance. We anticipate that these effects will be observable, and that 
understanding the physical mechanism will be essential to predict the behaviour of the composite materials. 

However a number of remarks related to our results are in order. First, the separation of high and low frequency 
regions depends on the natural scale qc oc 7/c (see Eqs. (^ - |l^)), and the large q 01 ui limit corresponds to the 
uncoupled case, when both eigenvectors are those of acoustic modes propagating more or less independently, while in 
the opposite limit the chains are tightly coupled. Of course in real materials the critical wave vector qc could be so 
large that the coupled limit is adequate over the whole accessible range of wave vectors. 

Second, in fact we considered a composite chain system embedded into two dimensional space. Thus one should 
distinguish between modes propagating along the x and y directions. For the modes propagating along x we found in 
the low q limit two diffusive and one pair of degenerate sound like modes. In addition there are transverse modes which 
are: two time-conjugated optical ones (i.e. with a gap) and two time-conjugated acoustic modes (sound like). For 
the propagation along y the situation is inverse: transverse modes give two diffusive modes and two time-conjugated 
sound like modes, while longitudinal modes are two time-conjugated optical modes and two time-conjugated acoustic 
modes. 

Note also that it is to be expected that the forces coupling the two chains are weak for relative displacements along 
the chain direction but of normal strength for displacements perpendicular to the chains. It implies that even the 
smallest wave vectors parallel to the x - axis probed e.g. by neutrons (see below our short survey of experimental 
data) may fall within the large q uncoupled regime, whereas all wave vectors in the orthogonal directions are in the 
normal coupled regime. 

Third, we assumed that atomic motions are not related to charge fluctuations (therefore we neglected Coulomb 
forces). As a result of Coulomb interactions additional q — gaps (plasmon gaps) will open, and in addition long-range 
Coulomb interactions will lead to a renormalization of the short-range elastic constants. The detailed discussion of 
all these problems is beyond the scope of our paper. 

It is worth to compare the predictions of the present model, with results of 0| where the authors also stressed the 
relevance of dynamical couplings. Qualitatively (especially in the low frequency region) our model leads to the same 
three longitudinal modes as predicted in 1^, 0| . However we presented a much simpler derivation which offers a deeper 
insight into the physical structure of the problem. Besides we should stress that we determined as well the dispersion 
laws for the transverse modes and for the high-frequency region we found an unusual g-independent contribution to 
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the sound modes attenuation (due to our dynamical coupling) superimposed on the conventional (oc q^) line- width. 
Note also that for our model there are no universal inequalities between sound mode velocities Ca,Cb,Cs as reported 
in 10]. 

In spite of the current research activity in the field of composite crystals, there are still few experimental results to 
which model predictions can be compared. Furthermore, the generic aspects of composite crystals which this and other 
models aim to describe, are often masked in practice by other effects, specific to each particular compound under study. 
One of the first composite systems studied with respect to dynamics, is the Hg^-gAsFg "mercury-chain" compound, 
which turned out to be a rather complex and unique system with two orthogonal sets of Hg - chains, showing Id 
liquid- like behaviour at high temperature {T > Tc = 120 A') Below Tc, the long-range phase-correlation between 
parallel Hg - chains is established through indirect interchain interactions, mediated by the orthogonal set of Hg - 
chains (which order simultaneously). Hence the compound is not a good Id - composite model system. 

The alkane-urea compounds appear a priori as better model systems because of their strongly uniaxial hexagonal 
structure. However, the structural transition at Tc ~ ISOiC, marking the onset of ferroelastic domains, precludes a 
detailed inelastic neutron or Brillouin study at low temperatures, where mode-broadening due to molecular disorder 
and anharmonic effects could be suppressed. The report by Schmicker et al., |jl^ of Brillouin scattering evidence for 
a propagating phase mode, has clearly been refuted by subsequent measurements by OUivier et al., ||20|| . Instead, 
these latter authors report observing essentially normal acoustic dispersions and a quasielastic line polarized in the 
chain direction, which they ascribe to an overdamped, pinned phase mode. In the limit when the mode damping 
coefficient / exceeds the pinning frequency ujg {ug ^ /), the phase- mode response would become quasielastic with 
a width jQ.E. = ^g/f- It is indeed quite plausible that, in the g-range covered by light scattering, neither ujg nor / 
would show much g-dependence, in which case one would not expect much g-dependence for the quasielastic linewidth 
iQ.Ei in agreement with the results in ref. ]20| . It conforms as well with our theoretical predictions presented in 
Section HI. The fact that 7q.e is also found to be T-independent between ISO-ftT and SOOA', is much more puzzling 
since, over such a wide T-range, one would expect ujg to soften and T to increase with T, leading to a significantly 
narrower quasielastic line at room-T, in clear contradiction with the measurements. Hence, the assignment of the 
observed quasielastic line as a collective excitation must be taken as tentative only. On the other hand inelastic 
neutron scattering data on the n=19 alkane-urea compound [2l[ |, show only one longitudinal acoustic (LA) mode 
in the chain direction (instead of two as would be expected in the high-q limit), with an anomalous q-independent 
linewidth. The non-observation of the collective modes associated with the alkane sublatticc can be ascribed to the 
strong orientational and translational disorder of the alkane molecules above Tc- 

Another interesting system is the aperiodic layered crystal Bi2Sr2CaCu20s+s, which was recently suggested to 
belong to the class of composite (rather than modulated) incommensurate systems [^5|. The mode structure 
observed by Etrillard et al. is consistent with the presence of two independent LA branches along the incommen- 
surate direction, as expected in the g-range of the neutron experiment. The relative strength of the two branches is 
however not understood, nor is the absence of a similar sublattice decoupling for the TA branches polarised in the 
incommensurate direction. 

It is obvious from the above remarks that the experimental situation still remains very open. To a large extent, this 
is due to the difhculty in identifying suitable model systems and in growing single-crystal specimens of appropriate 
size and quality for light and neutron scattering experiments. We hope that this work stimulates further efforts in 
this field. 

One of the authors (E.K) is indebted to RFFR and INTAS Grants for parial support. 

APPENDIX 

The equilibrium positions of the atoms in the incommensurate direction (x) could be written in a more abridged 
form than (^, namely 

x'^{t) = na + ul+ gi{na + - u^) ; 

xlXt) =mb + ul+ g2{mb + - <) , 

with gi {x + b) = gi (x) and 172 (x + a) = g2 (x) . Both functions gi and 172 can be found from a static microscopic model. 

The energy of interaction within the chains is written as usually, the dynamical coupling in the dissipation function 
can be written as: 

i? = ^(< - u„J^f{ma - nb) , 

mn 
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where f{x) can be modelled as: 



fix) = exp 



X 
"A2 



Suppose we have a relatively long-range coupling A ^ a, 6. In the case similar to that considered in the main text of 
the paper the equations of motion (after Fourier transformation over t) read 

-w V< = A„K+i + <_i - 2ul)+ 
i^^f{ma-nb)u''^; 

m 

+ir^/(ma-n6)u^, 

m 

or, after space Fourier transformation in ma, nb: = ?i°Q cxp ^(cjf — qna), = u^Qex.pi{ojt — qnb) 



■ 2 n 



2 6 



\ ■ b 



^ — (ma— n6)^/A^ Aqima—nh). a 



These equations give Xa^(^> ™) which we need for the analysis of the neutron scattering. 
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